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Background: Human papillomavirus (HPV)-positive oropharyngeal cancer (OPSCC) is associated with improved survival compared 
with HPV-negative disease. However, a minority of HPV-positive patients have poor prognosis. Currently, there is no generally 
accepted strategy for identifying these patients. 

Methods: We retrospectively analysed 270 consecutively treated OPSCC patients from three centres for effects of clinical, 
pathological, immunological, and molecular features on disease mortality. We used Cox regression to examine associations 
between factors and OPSCC death, and developed a prognostic model for 3-year mortality using logistic regression analysis. 

Results: Patients with HPV-positive tumours showed improved survival (hazard ratio (HR), 0.33 (0.21-0.53)). High levels of tumour- 
infiltrating lymphocytes (TILs) stratified HPV-positive patients into high-risk and low-risk groups (3-year survival; HPV-positive/ 
Tll_hj g h = 96%, HPV-positive/TIL| ow = 59%). Survival of HPV-positive/TIL| ow patients did not differ from HPV-negative patients (HR, 
1.01; P = 0.98). We developed a prognostic model for HPV-positive tumours using a 'training' cohort from one centre; the 
combination of TIL levels, heavy smoking, and T-stage were significant (AUROC = 0 • 87). This model was validated on patients 
from the other centres (detection rate 67%; false-positive rate 5.6%; AUROC = 0-82). 

Interpretation: Our data suggest that an immune response, reflected by TIL levels in the primary tumour, has an important role in 
the improved survival seen in most HPV-positive patients, and is relevant for the clinical evaluation of HPV-positive OPSCC. 
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Human papillomavirus (HPV) is the primary cause of oropha- 
ryngeal cancer (OPSCC), accounting for 40-80% of cases in the 
western world (D'Souza et al, 2007; Termine et al, 2008; Adelstein 
et al, 2009; Dayyani et al, 2010). The incidence of HPV-associated 
OPSCC has increased significantly over the past 40 years, and 
continues to rise (Adelstein et al, 2009; Chaturvedi et al, 2011). 

HPV-associated OPSCC commonly presents at advanced stage 
and is usually poorly differentiated (Adelstein et al, 2009; Sethi 
et al, 2012). Although both factors would be expected to result in 
poor outcomes, studies report significantly better long-term 
survival in most HPV-positive OPSCC patients compared with 
those with HPV-negative disease of matched stage (Licitra et al, 
2006; Fakhry et al, 2008; Worden et al, 2008; Lassen et al, 2009; 
Ang et al, 2010; Dayyani et al, 2010). The explanation for this 
remains unclear and may be multifactorial. Human papilloma- 
virus-positive tumours have improved response to chemotherapy/ 
radiotherapy compared with HPV-negative tumours, possibly 
modulated through the expression of wild- type p53 and Rb 
tumour suppressor genes; patients with HPV-positive cancers have 
a lower risk of developing secondary primary tumours, probably 
due to the absence of field cancerisation; in addition, HPV status 
inversely correlates with biomarkers of poor prognosis such as 
EGFR (Westra et al, 2008; Adelstein et al, 2009; Hong et al, 2010a). 

A potential explanation for survival benefit may be that tumour- 
infiltrating T lymphocytes (TILs), which permeate many HPV-positive 
OPSCC, confer a protective effect through an adaptive host immune 
response directed against viral antigens, and HPV- 16 E7-specific T 
cells have been reported in OPSCC patients (Albers et al, 2005; 
Wansom et al, 2010, 2012). Although TIL abundance is associated 
with improved clinical outcomes in a number of other tumour types, 
their role in head and neck malignancy is not well documented (Galon 
et al, 2006; Uppaluri et al, 2008; Gooden et al, 2011). 

Several studies have reported that the survival advantage of 
HPV-driven tumours is independent of treatment modality, 
whether surgery, radiation, or chemoradiotherapy (or combina- 
tions thereof), raising the possibility that some HPV-positive 
patients receive unnecessary treatment (Licitra et al, 2006; Fakhry 
et al, 2008; Adelstein et al, 2009; Lassen et al, 2009; Ang et al, 2010; 
Hong et al, 2010b). This has led to suggestions that treatment be 
de-intensified to reduce therapy- related morbidity (Adelstein et al, 
2009), and a number of clinical trials are testing this currently. 
However, within the HPV-positive OPSCC population, there 
remains a significant minority that respond poorly to treatment 
and have a poor prognosis. Ang et al (2010) identified three 
different 'at-risk' groups for patients with oropharyngeal cancer 
and found that HPV-positive patients who were heavy smokers 
and had advanced nodal disease (>10 pack-years smoking with 
N2b disease or above) were at an 'intermediate' risk of death. These 
factors have subsequently been validated independently and have 
been used as exclusion criteria in some de-escalation trials (Ang 
et al, 2010; Granata et al, 2012). Other authors have also described 
an association between heavy smoking and poor outcome in 
patients with HPV-positive tumours (Adelstein et al, 2009; Gillison 
et al, 2012; Sethi et al, 2012); however, there is still no widely 
accepted strategy for identifying high-risk HPV-positive patients. 

In this study, we assess the relative contributions of HPV status and 
pathological, molecular, and immunological features in survival from 
OPSCC. We correlate them with known clinical predictors of outcome 
to formulate a predictive tool for OPSCC patients with the ultimate 
aim of facilitating stratification of patients for treatment decisions. 



PATIENTS AND METHODS 



We retrospectively identified 442 consecutively treated patients 
with OPSCC (University Hospital Southampton (UHS; 2000-2010), 



n = 316; Poole NHS Foundation Trust and Bart's and the 
London NHS Trust (PFT/BLT; 2000-2006), n = 126. Patients were 
excluded from further analysis if they had follow-up of <6 
months, had received either no treatment or treatment without 
curative intent, or had an unknown cause of death. Following these 
exclusions, 402 patients remained. Archived pathology material 
was available for 274 patients who were subsequently analysed 
further. Demographics and survival of patients without pathology 
material were comparable to those with pathology material 
(Supplementary Information). All patients were analysed accord- 
ing to REMARK guidelines (McShane et al, 2005). Data were 
obtained from patient records on age, sex, site, smoking (non-, 
ex- (stopped smoking >6 months previously), light smokers 
(defined as < 10 pack year history (Ang et al, 2010; Gillison et al, 
2012)), heavy smokers ( > 10 pack year history)), alcohol (non-/ex- 
drinker, current drinker > 10 units per week), T stage, N stage, 
treatment modality, treatment failure, tumour recurrence, date, 
and cause of death. Data on smoking and alcohol consumption can 
be difficult to accurately obtain in retrospective studies, often 
because of poor documentation. In this study, smoking data were 
more accurately recorded in the medical records than alcohol 
consumption, with reliable data available in 87% and 74% of cases, 
respectively. Demographics are shown in Table 1. Ethical approval 
was obtained (UKCRN 8130; ISRCTN 71276356; and REC 
references 09/H0501/90 and 07/Q0405/1). 

Treatment. A treatment summary is included in Table 1. The 
centres involved in this study differed in their approach to 
treatment of OPSCC, with more patients being surgically treated at 
UHS. This allowed a direct comparison of survival according to 
treatment modality. Patients undergoing chemoradiotherapy 
received 2-6 doses of platinum-based chemotherapy (either 
cisplatin or carboplatin) combined with 64-66 Gy of radiotherapy 
given in 32-33 fractions. No patients received accelerated 
fractionation or intensity-modulated radiotherapy. Neoadjuvant 
chemotherapy was 1-2 doses of either cisplatin or carboplatin, in 
conjunction with 5-flurouracil. Of the surgically treated patients, 
78% received post- operative radiotherapy. Full details are given in 
Supplementary Information. 

Histology and immunochemistry. Tumour histology was 
reviewed by pathologists (GJT, SJ, and KP) blinded to treatment 
and outcome, according to the 1998 UK Royal College of 
Pathologist Guidelines, and a suitable paraffin block was selected. 
From a haematoxylin/eosin (H + E) -stained section, data were 
obtained on grade (well, moderate, or poorly differentiated). 
A prominent lymphocytic infiltrate (TIL) was scored under 
low-power magnification ( x 2 • 5 objective) as high (diffuse; 
present in >80% of tumour/stroma), low (weak/absent; present 
in <20% of tumour/stroma), or moderate (patchy; present in 
20-80% of tumour/stroma) as described (Marsh et al, 2011). 
Concordance between pathologists in identifying TILi ow tumours 
was > 95%. Typical examples of TIL high and TIL low tumours are 
shown in Figure ID and Supplementary Information. 

Tissue microarrays (TMAs) were constructed using triplicate, 
randomly selected, paraffin -embedded 1-mm tumour cores 
(Aphelys Minicore 2, Mitogen, Harpenden, UK). Automated 
immunostaining (Ventana XT, Ventana, Tucson, AZ, USA) was 
performed in a CPA-accredited clinical cellular pathology depart- 
ment using antibodies optimised to national diagnostic standards 
(NEQAS). The antigenic targets were as follows: pl6 (CINtec, 
Roche, Burgess Hill, UK), EGFR (Dako, Ely, UK), T-lymphocyte 
markers (CD3, pan T cells; CD4, helper T cells; CD8, cytotoxic T 
cells (all from Novocastra, Milton Keynes, UK); FoxP3, regulatory 
T cells (eBioscience, Hatfield, UK), and B lymphocytes (CD20) 
(Novocastra). TILs were quantified using a Zeiss AxioCam MRc5 
microscope (Zeiss, Cambridge, UK) and Zeiss Axiovision software 
(version 4 • 8 • 1 • 0) (Zeiss). An average intratumoural TIL 
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Table 1. Baseline characteristics of patients in 1 


he retrospective patient cohorts (i.e., at the time of diagnosis) 








All OPSCC 

n = 274 
Number (%) 


HPV-positive OPSCC 
n = 149 
Number (%) 


HPV-negative OPSCC 
n=121 
Number (%) 


Gender 


Male 
Female 


200 (73.0) 
74 (27.0) 


111 (74.5) 
38 (25.5) 


86 (71.1) 
35 (28.9) 


Age at diagnosis (years) 


<50 
50-69 
70 + 

Mean (s.d.) 


58 (21.2) 
168 (61.3) 

48 (17.5) 
58.2 (11.2) 


40 (26.8) 
91 (61.1) 
18 (12.1) 
56.4 (10.8) 


18 (14.9) 
74 (61.2) 
29 (24.0) 
60.4 (11.4) 


Smoking 


Non-/ex-smoker 
Current <10 pack-years 
Current >10 pack-years 
Not Known 


107 (39.1) 
18 (6.6) 

114 (41.6) 
35 (12.8) 


82 (55.0) 
10 (6.7) 
41 (27.5) 
16 (10.7) 


23 (19.0) 

8 (6.6) 
72 (59.5) 
18 (14.9) 


Alcohol 


Non-/ex-drinker 
Current drinker 
Not known 


34 (12.4) 
167 (60.9) 
73 (26.6) 


20 (13.4) 
94 (63.1) 
35 (23.5) 


14 (11.6) 
72 (59.5) 
35 (28.9) 


Tumour site 


Tonsil 

Tongue base 
Other oropharynx 


158 (57.7) 
70 (25.5) 
46 (16.8) 


99 (66.4) 
40 (26.8) 
10 (6.7) 


57 (47.1) 
28 (23.1) 
36 (29.8) 


Median length of follow-up in years (range) 


58.0 (8-137) 


58.5 (8-137) 


57 (8-128) 


Disease stage 


l/ll 
lll/IV 

Not known 


53 (19.3) 
219 (79.9) 
2 (0.7) 


10 (6.7) 
138 (92.6) 
1 (0.7) 


42 (34.7) 
78 (64.5) 
1 (0.8) 


T stage 


T1/2 
T3/4 

Not known 


187 (68.2) 
82 (29.9) 
5(1.8) 


107 (71.8) 
41 (27.5) 
1 (0.7) 


77 (63.6) 
41 (33.9) 
3 (2.5) 


Nodal metastases 


No 
Yes 

Not known 


68 (24.8) 
204 (74.5) 
2 (0.7) 


15 (10.1) 
133 (89.3) 
1 (0.7) 


52 (43.0) 
68 (56.2) 
1 (0.8) 


N stage 


N0-N2a 
N2b-N3 
Not known 


142 (51.8) 
130 (47.4) 
2 (0.7) 


68 (45.6) 
80 (53.7) 
1 (0.7) 


73 (60.3) 
47 (38.8) 
1 (0.8) 


Grade 


Well/moderately differentiated 
Poorly differentiated 


103 (37.6) 
171 (62.4) 


34 (22.8) 
115 (77.2) 


68 (56.2) 
53 (43.8) 


Treatment 


Surgery ± radiotherapy 

Radiotherapy 

Chemoradiotherapy 


117 (42.7) 
59 (21.5) 
98 (35.8) 


60 (40.3) 
24 (16.1) 
65 (43.6) 


55 (45.5) 
34 (28.1) 
32 (26.4) 


TIL status 


Low 

Moderate 
High 

Not known 


79 (28.8) 
101 (36.9) 
92 (33.6) 
2 (0.7) 


22 (14.8) 
53 (35.6) 
72 (48.3) 
2(1.3) 


56 (46.3) 
46 (38.0) 
19 (15.7) 
0 (0.0) 
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Table 1. (Contin 



All OPSCC 

n = 274 
Number (%) 



HPV-positive OPSCC 
n = 149 
Number (%) 



HPV-negative OPSCC 
n=121 
Number (%) 



HPV status 

Negative 
Positive 
Not known 



121 (44.2) 
149 (54.4) 
4(1.5) 



EGFR 



Negative 
Positive 



108 (39.4) 
166 (60.6) 



77 (51.7) 
72 (48.3) 



29 (24.0) 
92 (76.0) 



Abbreviations: HPV = human papillomavirus; OPSCC = oropharyngeal cancer; TIL = tumour-infiltrating lymphocyte. 



score/high-power field ( x 400) was calculated across representa- 
tive areas of each triplicate core to allow for intratumoural 
heterogeneity (i.e., an average of three counts across the tumour). 
EGFR and pl6 were scored as positive (>50% tumour cells 
positive) or negative (^50% tumour cells positive). 

In situ hybridisation and HPV classification. HPV in situ 
hybridisation (ISH) was performed using Genpoint HPV DNA 
Probe Cocktail (Dako), which recognises the 13 most prevalent 
high-risk HPV genotypes, including HPV- 16 and -18. Cases were 
scored as negative or positive. Tumours positive both for pl6 
immunochemistry and HPV ISH were classified as HPV positive 
(Robinson et al, 2010). 

Microarray data analysis. Affymetrix U133 Plus 2.0 CEL files for 
46 HNSCC patient samples (24 survivors and 22 non-survivors) 
were downloaded from a publicly available microarray data set 
(http://bioinformatics.picr.man.ac.uk/vice/PublicProject.vice?pid=361) 
and normalised in BioConductor-R (version 2.15.1) (Bioconductor, 
Seattle, WA, USA). Gene-set enrichment analysis was performed 
using the GSEA 2.7 software package (Broad Institute of MIT and 
Harvard, Cambridge, MA, USA). All 15,031 genes differentially 
expressed between survivors and non- survivors were rank ordered 
and tested for enrichment of gene sets upregulated during T-cell 
differentiation (Best et al, 2013). 

Statistics. The primary endpoint was death from OPSCC (i.e., 
disease- specific survival (DSS)). An OPSCC-specific survival time 
was measured from the date of diagnosis until date of death from 
OPSCC or date last seen alive. For those that died of other causes, 
survival from OPSCC was censored at the time of death. Patients in 
whom the cause of death was unknown, or who received either no 
treatment or treatment with palliative intent, were excluded from 
all survival analyses. Kaplan-Meier plots (with log-rank tests) and 
both unadjusted and adjusted Cox proportional hazards models 
were used to analyse the survival data. Correlations between factors 
were assessed by cross-tabulation and the use of either Fisher's 
exact test or # 2 -test where appropriate. 

To evaluate the prognostic performance of each factor (i.e., the 
3 -year OPSCC death rate), we calculated the detection rate (DR; 
also known as sensitivity) and false-positive rate (FPR; also known 
as 1 — specificity). Detection rate is the proportion of patients who 
died from OPSCC with marker-positive results; FPR is the 
proportion of patients who did not die from OPSCC with 
marker-positive results. Likelihood ratios (DR/FPR) were obtained, 
which indicate the strength of a marker (i.e., maximising the DR 
while minimising the FPR; Deeks and Altman, 2004). 

A prognostic model for 3 -year mortality was developed using 
binary logistic regression, with a backwards selection procedure 
(Royston et al, 2009). Variables were retained in the model if they 



were significant at the 5% level. The UHS cohort was used as the 
'training' set to develop the model, with the PFT/BLT cohort acting 
as the Validation' set. A Receiver- operating characteristic (ROC) 
curve of the linear predictor was used to identify the optimum cut 
point (balancing both sensitivity and specificity; Royston et al, 
2009). The prognostic ability of the final model was then assessed 
through the calculation of DR and FPR values in the Validation' set 
(Heagerty et al, 2000). All statistical tests were two sided and a 
P- value of 0.05 was considered significant. Statistical analyses 
were performed using SPSS for Windows version 19.0 (IBM, 
Portsmouth, UK). 



RESULTS 



Clinical, pathological, and immunochemical characteristics at the 
time of diagnosis are shown in Table 1. Human papillomavirus - 
positive tumours were defined as those positive for both pl6 and 
HPV ISH (Robinson et al, 2010). Similar to previous studies, 94% 
of pl6-positive cases were HPV ISH positive and 89% of HPV ISH- 
positive cases were pl6 positive (Robinson et al, 2010). 

Of the 274 patients included in the study, 149 (54.4%) were 
classified as HPV positive. Human papillomavirus status could not 
be established in four cases (insufficient tumour tissue remaining 
in the paraffin block); thus, further analysis was performed on 270 
patients (median follow-up 58 months and minimum follow-up 8 
months). There were 82 deaths from OPSCC over the study period. 
Consistent with previous studies, patients with HPV-positive 
tumours had significantly improved DSS compared with those 
that were HPV negative (Figure 1 A; hazard ratio (HR) = 0.33 
(0.21-0.53)). Table 2A shows the HRs for OPSCC mortality for 
each factor separately. Improved survival was significantly 
associated with HPV-positive tumours, a younger age at diagnosis, 
non-, ex-, or light smokers (defined as < 10 pack year history (Ang 
et al, 2010; Gillison et al, 2012)), early T-stage tumours (Tl/2), 
early N-stage tumours (N0-N2a), tumours that were EGFR 
negative, and a high level of TILs (assessed on an H + E-stained 
section). When a multivariate analysis was performed (adjusting 
for age, T stage, N stage, and smoking status), HPV status and TIL 
levels remained statistically significant markers of OPSCC 
mortality (Table 2A). 

In unstratified OPSCC, high TIL levels predicted for survival 
(Figure 1C; Table 2A) and correlated significantly with HPV- 
positive tumours (P< 0.001). The percentage of HPV-positive 
patients with high, moderate, and low TIL was 49%, 36%, and 
15% respectively (16%, 38%, 46%, respectively, in HPV- 
negative tumours). Figure ID and F show typical examples of 
HPV-positive tumours with high and low TIL levels. 
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No. at risk 
HPV + /TIL high 
HPV+/TIL mod 53 38 
HPV+/TIL low 22 12 
HPV- 



Years 



72 59 39 16 



28 



11 

4 



10 12 



121 67 36 20 




































pi 6 


x100 


CD3 






Figure 1 . Kaplan-Meier curves for OPSCC mortality in combined patient cohorts. (A) Kaplan-Meier curves for OPSCC mortality according to HPV 
status (log-rank test, P< 0-001). (B) Immunochemistry (p16) and ISH (HPV) showing an example of an HPV-positive tumour. (C) Kaplan-Meier 
curves for OPSCC mortality according to abundance of infiltrating lymphocytes (TIL; high, moderate or low; log-rank test, P< 0-001). (D) H&E- 
stained sections showing examples of OPSCC with high and low TIL infiltrates. Tumour islands are arrowed. In TIL^h tumours, lymphocytes fill the 
stromal compartment and infiltrate into tumour islands. (E) Kaplan-Meier curves for OPSCC mortality stratified according to HPV status and TIL 
levels. Patients with HPV + /TIL| ow tumours show similar survival to HPV-negative patients. (F) Immunochemistry showing typical examples of HPV- 
positive OPSCC with high and low TIL. Tumour is positive for p16, T lymphocytes for CD3. 



Demographics of these groups are shown in Supplementary 
Information. 

There was a highly significant difference between the 
Kaplan-Meier curves for HPV-positive OPSCC DSS according to 
TIL levels (P< 0.001; Figure IE); patients with HPV-positive/TILi ow 
tumours showed similar DSS to HPV-negative patients (Figure IE) 
with an adjusted HR of 1.01 (P = 0.98; Table 2A). TIL levels 
predicted similarly for progression-free survival (PFS; log-rank 
P< 0.001; data not shown); patients with HPV-positive/TIL low 
tumours also showed similar PFS to HPV-negative patients with an 
adjusted HR of 0.64 (P = 0.23). The adjusted HRs for HPV-positive 



tumours with moderate and high TIL levels indicated a 49% and 
79% reduction in OPSCC mortality compared with HPV-negative 
tumours (adjusted HRs, 0-51 and 0-21, respectively; P = 0.17 and 
P = 0.02; Table 2B). In HPV-negative tumours, patients with TIL high 
tumours did not show significantly improved survival (adjusted HR, 
0.47, P = 0.17) over those with TIL low tumours. 

Similar to previous studies, we found that heavy smoking 
reduced survival benefit in HPV-positive patients (HR, 5.81 
(2.38-14.16)). We therefore assessed whether there was an 
association between heavy smoking and low TIL levels; low TIL 
levels did not correlate with heavy smoking ( > 10 pack-years) and 
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Table 2A. HRs for disease-specific death from OPSCC 






All OPSCC (n = 274) 




Univariate 
HR (95% CI) 


P- value 


Multivariate HR a 
(95% CI) 


P- value 


Age 


For each additional year 


1 .02 (1 .00-1 .04) 


0.04 


1.03 (1.01-1.05) 


0.01 


Gender 


Female 
Male 


1 

1.49 (0.89-2.48) 


0.13 


1 

1.10 (0.61-1.96) 


0.75 


Smoking 


Non/ex-smoker/current smoker <10 pack-years 
Current smoker >10 pack/years 


1 

3.24 (1.94-5.42) 


<0.001 


1 

3.38 (2.01-5.70) 


<0.001 


Drinking 


N on/ex-drinker 
Current drinker 


1 

0.67 (0.36-1.24) 


0.20 


1 

0.44 (0.22-0.88) 


0.02 


Stage 


Early (l/ll) 
Late (I I l/l V) 


1 

1.01 (0.58-1.74) 


0.98 






T stage 


T1/2 
T3/4 


1 

1.89 (1.21-2.94) 


0.005 


1 

1.56 (0.95-2.57) 


0.08 


Nodal metastases 


No 
Yes 


1 

1.11 (0.67-1.84) 


0.68 






N stage 


N0-N2a 
N2b-N3 


1 

1.68 (1.08-2.60) 


0.02 


1 

1.75 (1.07-2.86) 


0.03 


Tumour grade 


Well/moderately differentiated 
Poorly differentiated 


1 

0.80 (0.52-1.24) 


0.31 


1 

0.97 (0.59-1 .60) 


0.90 


Treatment 


Surgery 

Radiotherapy 

Chemoradiotherapy 


1 

1.53 (0.90-2.61) 
0.95 (0.57-1.58) 


0.12 
0.83 


1 

1.06 (0.56-2.03) 
0.91 (0.51-1.63) 


0.85 
0.74 


HPV status 


Negative 
Positive 


1 

0.33 (0.21-0.53) 


<0.001 


1 

0.45 (0.26-0.79) 


0.005 


TIL Status 


TIL low 

TIL moderate 

TIL high 


1 

0.59 (0.37-0.95) 
0.21 (0.11-0.41) 


0.03 
<0.001 


1 

0.66 (0.39-1.13) 
0.28 (0.13-0.62) 


0.13 
0.002 


HPV/TIL combined 


HPV- 

HPV + /TIL low 
HPV + /TIL moderate 
HPV + /TIL high 


1 

0.83 (0.40-1.75) 
0.41 (0.22-0.78) 
0.14 (0.06-0.33) 


0.63 
0.006 
<0.001 


1 

1.01 (0.45-2.24) 
0.49 (0.24-1.01) 
0.17 (0.06-0.49) 


0.98 

0.054 

0.001 


Abbreviations: CI = confidence interval; HPV= human papillomavirus; HR = hazards ratio; OPSCC = oropharyngeal cancer; TIL = tumour-infiltrating lymphocyte. 
Multivariate analysis adjusted for age, T stage, N stage, and smoking status. Patients with no treatment excluded from analysis. 
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there was no significant difference in the proportion of TIL low 
tumours according to the amount smoked (percentage of TIL low 
tumours: <10 pack year history 14.1%, > 10 pack year history 
22.5%, P = 0.15). 

We performed more complex analysis of the immune infiltrate 
to establish whether this improved prognostication in HPV- 
positive tumours. We quantified intratumoural T-lymphocyte 
number (CD3) and T-cell subsets (CD4, CD8, and Foxp3). We 
performed ROC analysis (death from OPSCC at 3 years) to 
determine the predictive value of T-cell number and/or subtype, 
including subtype ratios (CD4:CD8 and CD8:FoxP3), in compar- 
ison with simple H + E evaluation of TIL infiltrate. Immuno- 
chemical quantification did not outperform the simpler scoring 
method (Figure 2 A and B, and Supplementary Information). 

To evaluate whether high TIL levels reflected a bystander or 
effector response, we tested a publicly available HNSCC microarray 
data set (Thurlow et al, 2010) for enrichment of gene sets 
defined as upregulated in naive, effector, or memory T cells 
(Best et al, 2013). High expression of gene sets associated with 
CD8 + T-cell effector function were significantly enriched 
(enrichment score (ES) = 0.70, P = 0.046) in patients with better 
survival (Figure 3 and Supplementary Information). Genes 
upregulated in naive or late memory CD8 + T cells were 
overrepresented in patients with poor survival but at non- 
significant levels, possibly reflecting their TIL low status 
(ES = 0.30, P = 0.46; data not shown). 

Table 3 shows the prognostic performance of each factor (using 
DR and FPR) for predicting 3-year OPSCC mortality in HPV- 
positive and HPV-negative patients. Commonly reported prog- 
nostic factors are included. In HPV-positive OPSCC, TIL levels 
were the best-performing factor, with a likelihood ratio (LR) for 
TIL low tumour of 3.3, followed by heavy smoking and advanced 
T stage (LR 2.78 and 2.44, respectively). In HPV-negative OPSCC, 
advancing age, advanced stage, and poor differentiation were the 
predictive factors with the highest LRs (Table 3). 

HPV-positive patients had received different types of therapy 
(cohorts combined: surgery n = 55, chemoradio therapy n = 65, 
and radiotherapy n = 24). When adjusted for stage and age, there 
were no significant treatment- dependent survival benefits (HR for 
DSS relative to surgery: radiotherapy 0-76, P = 0.79 and 
chemoradiotherapy 1 • 25, P= 0.84). 

Using the UHS cohort as a training set, we used logistic 
regression to develop a predictive model for survival from HPV- 
positive OPSCC (at 3 years). The initial variables entered into the 
model were as follows: age, T stage, N stage, grade, HPV status, 
TIL level, alcohol history, EGFR, and smoking history. The final 
model included TIL levels (low vs moderate or high), T stage (Tl/2 
vs T3/4), and smoking status (heavy vs light-, ex-, or non-smoker). 
The regression parameter estimates were used to create a 
prognostic score (equation: — 4.948 + 2.768 (if a current heavy 
smoker) + 2.310 (if TILi ow ) + 2.928 (if T3/4)). ROC analysis was 
performed and produced an area under the curve of 0-87 
(Figure 4A). A cut-off score of — 0.945 was selected based on 
this ROC curve (DR, 72.7%; FPR 10.2%; high risk, n = U) and 
applied to the PFT/BLT (test; high risk, n = 7) cohort. This score 
was highly predictive of 3 -year mortality, with 66.7% of dead 
patients test-positive (DR = 67%) and 5.6% of alive patients test- 
positive (FPR = 5.6%). The resultant LR for prediction of 3-year 
mortality was 11.9 (PFT/BLT, AUROC, 0-82; Figure 4B). Repla- 
cing TIL levels with counts or ratios of T-lymphocyte subsets did 
not improve this model. 



DISCUSSION 



Patients with HPV-positive OPSCC have significantly better 
survival, although the biological reason for this remains unclear 



(Adelstein et al, 2009). To reduce therapy-associated morbidity, it 
may be possible to de-intensify treatment, particularly as several 
studies have shown survival benefit to be independent of therapy 
choice (Licitra et al, 2006; Fakhry et al, 2008; Worden et al, 2008; 
Adelstein et al, 2009; Lassen et al, 2009; Ang et al, 2010; Dayyani 
et al, 2010; Hong et al, 2010b). However, a significant minority of 
HPV-positive OPSCC patients experience treatment failure and 
have poor outcome; it has been suggested that future clinical 
investigations concentrate on this group in order to improve 
outcome (and also to avoid treatment de-intensification in these 
patients; Adelstein et al, 2009). Currently, there is no widely 
accepted strategy for the identification of poor prognosis in HPV- 
positive patients before therapy. 

A plausible explanation for the differences in survival between 
HPV-positive and HPV-negative patients is that virally 
driven tumours provoke an adaptive immune response directed 
against tumour- expressed viral antigens: immune responses against 
the foreign viral antigen are less likely to be suppressed; there is no 
central immunological tolerance to confound the immune system's 
attempts to control cancer. HPV-16-specific CD8-positive 
T cells have been detected in the blood of HPV-positive OPSCC 
patients and, more recently, isolated from tumours, implicating a 
role in the anti-tumour response (Albers et al, 2005; Heusinkveld 
et al, 2012; Wansom et al, 2012). Several recent studies have 
further demonstrated anti-tumour immunity in HPV-positive 
OPSCC. Liang et al (2012) found better survival in HPV-positive 
patients with seropositivity to E6/E7. Furthermore, infiltration 
of HPV-positive HNSCC by PD-1 -expressing T lymphocytes 
has been shown to be a favourable prognostic factor 
(Badoual et al, 2013). 

In several solid tumours, quantification of densities of various 
T-cell subsets, including CD8 + , CD4 + Treg cells (FoxP3 + 
CD4 + ), or subset ratios, have been suggested to improve predictive 
power over that of absolute T-cell number (CD3 + ; Gooden et al, 
2011). Cytotoxic CD8 + T cells are the principal anti-tumour 
effector cells, and their abundance has been shown to be a 
predictor of positive outcome in several tumour types, particularly 
colorectal cancer, suggesting that the adaptive immune system has 
a role in suppressing tumour progression (Galon et al, 2006). Our 
analysis of a publicly available HNSCC microarray data set showed 
enrichment of genes associated with CD8 + T-cell effector function 
in patients with improved survival, arguing that the T cells are not 
simply 'innocent bystanders' of the biological events in the tumour 
but take an active role in tumour recognition. 

We conducted this study to examine the effects of HPV status 
and TIL levels on survival in OPSCC patients. Similar to other 
large studies, HPV-positive tumours were associated with sig- 
nificantly improved survival (3-year survival; HPV-positive 82% vs 
HPV-negative 56%, P< 0.001); pathological and clinical features of 
HPV-positive tumours were consistent with these series. The HPV- 
positive OPSCC were generally of late stage and were poorly 
differentiated (both P< 0.001). HPV-positivity inversely correlated 
with smoking, and with EGFR expression (all P< 0.005). TIL levels 
also correlated significantly with the HPV status (P< 0.001), with 
around 85% of HPV-positive OPSCC tumours containing high or 
moderate levels of TIL. Taking into account the array data, this is 
likely to reflect an adaptive anti-tumour response as the 
mechanism for improved outcome. TIL levels stratified HPV- 
positive OPSCC patients into those with good and poor prognosis. 
The 3 -year survival for HPV-positive/TIL high tumours was 96% 
compared with 76% for HPV-positive/TIL mod and 59% for HPV- 
positive/TIL low . The HPV-positive/TIL low tumours had similar 
survival to HPV-negative tumours (3-year survival, 56%). We also 
quantified densities of CD3-, CD4-, CD8-, and FoxP3 -positive T 
cells and carried out ROC analysis to determine the predictive 
value of cell numbers or subset ratios. These more complex 
analyses did not outperform the simpler scoring method 
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Table 2B. HRs for disease-specific death in OPSCC stratified by HPV status 



HPV-positive OPSCC # n = 1 49 " HPV-negative OPSCC, n = 1 21 





Univariate HR 
(95% CI) 


Multivariate 
HR a 
(95% CI) 


P- value 


Univariate 

HR 
(95% CI) 


Multivariate 
HR a 
(95% CI) 


P- value 


Age 


For each additional year 


1 .00 (0.97-1 .04) 


1 .02 (0.98-1 .06) 


0.28 


1 .02 (0.99-1 .04) 


1 .02 (0.99-1 .05) 


0.19 


Gender 


Female 
Male 


1 

1.19 (0.50-2.81) 


1 

0.74 (0.26-2.11) 


0.58 


1 

2.06 (1 .06-4.02) 


1 

1.70 (0.79-3.64) 


0.17 


Smoking 


N on/ex-smoker/current Smoker <10 pack-years 
Current smoker >10 pack-years 


1 

5.81 (2.38-14.16) 


1 

5.80 (2.31-14.54) 


<0.001 


1 

1.37 (0.71-2.64) 


1 

1.61 (0.82-3.16) 


0.17 


Drinking 


N on/ex-drinker 
Current drinker 


1 

0.63 (0.23-1.71) 


1 

0.29 (0.08-1 .04) 


0.06 


1 

0.68 (0.31-1.49) 


1 

0.54 (0.22-1.28) 


0.16 


Stage 


Early (l/ll) 
Late (I I l/l V) 


1 

1.93 (0.26-14.27) 






1 

1.83 (0.99-3.37) 






T stage 


T1/2 
T3/4 


1 

3.42 (1.61-7.28) 


1 

3.04 (1.28-7.22) 


0.01 


1 

1.24 (0.71-2.17) 


1 

0.98 (0.52-1 .84) 


0.94 


Nodal metastases 


No 

Yes 


1 

1.65 (0.39-7.00) 






1 

2.13 (1.19-3.79) 






N stage 


N0-N2a 
N2b-N3 


1 

2.29 (1.00-5.24) 


1 

1.67 (0.67-4.14) 


0.27 


1 

2.06 (1.20-3.53) 


1 

2.15 (1.14-4.06) 


0.02 


Tumour grade 


Well/moderately differentiated 
Poorly differentiated 


1 

0.74 (0.32-1.69) 


1 

0.74 (0.30-1 .84) 


0.52 


1 

1.58 (0.92-2.72) 


1 

1.50 (0.81-2.81) 


0.20 


Treatment 














Surgery 

Radiotherapy 

Chemoradiotherapy 


1 

1.51 (0.55-4.15) 
0.98 (0.42-2.31) 


1 

1.01 (0.30-3.36) 

1 .02 (0.39-2.70) 


0.99 
0.96 


1 

1.64 (0.87-3.12) 
1.48 (0.76-2.88) 


1 

1.50 (0.70-3.22) 
1.18 (0.55-2.57) 


0.30 
0.67 


TIL status 


TIL high 

TIL moderate 

TIL low 


1 

2.96 (1.11-7.89) 
5.67 (1.96-16.38) 


1 

2.47 (0.72-8.50) 
4.86 (1.34-17.60) 


0.15 
0.02 


1 

1.39 (0.55-3.48) 
1.71 (0.71-4.13) 


1 

1.59 (0.54-4.71) 
2.11 (0.73-6.16) 


0.40 
0.17 


EGFR 


Negative 
Positive 


1 

1.65 (0.76-3.62) 


1 

1.54 (0.63-3.74) 


0.34 


1 

0.99 (0.53-1.87) 


1 

1.31 (0.61-2.83) 


0.49 


Abbreviations: CI = confidence interval; HPV=human papillomavirus; HR = hazards ratio; OPSCC = oropharyngeal cancer; TIL = 
Multivariate analysis adjusted for age, T stage, N stage, and smoking status. Patients with no treatment excluded from analys 


= tumour-infiltrating lymphocyte, 
s. 





performed on H + E-stained sections. It is possible that lympho- 
cyte counts performed on randomly selected TMA cores may not 
be representative of the tumour as a whole. However, TMA cell 
counts significantly correlated with whole-section TIL levels, 



arguing against this (e.g., total (CD3 + ) T-lymphocyte count: 
TILi ow 42.9, TIL mod 60.4, and TIL high 96.3, P< 0.001). 

In keeping with previously published studies (Badoual et al> 
2006; Bron et al, 2013), we found that infiltration with 
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TIL levels 


0.76 


Total T-cell count 


0.70 
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0.64 
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0.70 
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CD4:CD8 ratio 


0.41 
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0 0.2 0.4 0.6 0.8 1.0 
False-positive rate (1 -specificity) 

Figure 2. Receiver-operating characteristic (ROC) curve for 3-year mortality in HPV-positive OPSCC to determine the predictive value of T-cell 
number and/or subtype. (A) Immunochemistry showing representative example of an OPSCC stained with T-cell markers (CD3, pan T cells; CD4, 
helper T cells; CD8 cytotoxic T cells; and FoxP3, regulatory T cells). (B) ROC curve for 3-year HPV-positive OPSCC mortality using TIL levels (scored 
on H + E-stained section), compared with cell counts following immunochemistry for T-cell subtypes (CD3 + , CD4 + , CD8 + , and Foxp3 + ), 
including subtype ratios (CD4:CD8, and Foxp3:CD8). Immunochemical quantification did not outperform the simpler scoring method. 



Enrichment plot: Effector. grp 
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Figure 3. T cell gene-set enrichment in HNSCC. Gene-set enrichment 
analysis of publicly available gene sets describing effector CD8 + T 
cells. Genes to the left and right of the rank-ordered list are enriched in 
HNSCC patients with early (survival <3 years, mean 16 months) or late 
mortality (survival >3 years, mean 49 months), respectively. The rank- 
ordered list of genes describing effector CD8+ T cells is shown in 
Supplementary Information. 

FoxP3 -positive regulatory T cells was associated with improved 
survival (Figure 2B). However, we feel that the increased numbers 
seen in HPV-positive tumours simply reflect an increase in overall 
T-lymphocyte numbers in these tumours. Indeed, the proportion 
of Foxp3 -positive cells was reduced in HPV-positive tumours, 
compared with those that were HPV-negative (Supplementary 



Information). Interestingly, TIL levels did not predict for survival 
in HPV-negative tumours. 

Several studies have reported potential clinical and molecular 
predictors for identifying poor prognosis in patients with HPV- 
positive tumours. The HPV-positive light- or non-smokers have 
the best outcome in OPSCC (Hafkamp et al, 2008; Kumar et al, 
2008; Gillison et al, 2012; Sethi et al, 2012). We also found that 
HPV-positive current heavy smokers had reduced survival benefit 
compared with non-, ex-, or light smokers (3-year survival; HPV- 
positive/current heavy smoker 63% vs HPV-positive/non-, 
ex-, current light smoker 94%; Supplementary Information). There 
was no significant association between heavy smoking and low TIL 
levels, arguing against TIL levels acting as a surrogate marker for 
heavy smoking. However, it is possible that poor documentation in 
retrospective records may be responsible for this lack of 
correlation, and this should be examined prospectively. We 
calculated LRs to indicate the strength of a marker (maximising 
the DR while minimising the FPR). In HPV-positive patients, low 
TIL levels and smoking status were the best predictors of 
mortality at 3 years (LRs 3-3 and 2-8, respectively). This was 
reflected in our prognostic model where, in the 'training' set 
(UHS patients) only TIL levels, smoking, and T stage were 
significant (AUROC analysis 0 • 87). When this model was 
applied to the Validation set' (PFT/BLT cohort), it was found 
to be highly predictive (DR 66.7%; FPR 5.6%; LR for 3-year 
mortality 11.9). This LR is better than other currently used 
predictive tests. For example, the unrehydrated guaiac fecal occult 
blood test used in screening for colorectal carcinoma has an LR of 
6.4 (DR 64% and FPR 10%) (Allison et al, 2007). The 3-year 
survival of low-risk HPV-positive patients (score < — 0.945) was 
94% compared with 36.8% for those that were classified as 
patients with high risk (score > - 0.945, P< 0.001). 

Ang et al (2010) have suggested a predictive model for HPV- 
positive OPSCC where patients who are heavy smokers ( > 10 pack 
year history) and have advanced nodal disease (>N2b) are at 
increased risk of death. We compared this model with our own, 
and calculated both LRs and 3 -year survival rates. The LR for the 
Ang model in our combined cohort was 4.8 (DR 61.1% and FPR 
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Table 3. LRs for 3-year death from OPSCC 






















Prognostic values for prediction of death from OPSCC at 3 years 






HPV-positive OPSCC 


HPV-negative OPSCC 




DR % (21 deaths) 


FPR % (95 alive) 


LR 


DR % (41 deaths) 


FPR % (53 alive) 


LR 


Age # years 


^50 


57.1 


72.6 


0.79 


85.4 


81.1 


1.05 


^55 


47.6 


51.6 


0.92 


73.2 


58.5 


1.25 


^60 


47.6 


30.5 


1.56 


61.0 


35.8 


1.70 


Gender 


Male 


81.0 


66.3 


1.22 


85.4 


62.3 


1.37 


Smoking 


Smoker >10/day 


77.8 


27.9 


2.78 


71.4 


66.7 


1.07 


Alcohol 


Drinker 


72.2 


86.1 


0.84 


75.9 


87.5 


0.86 


Stage 


^W, III, IV 


100.0 


97.9 


1.02 


97.6 


73.6 


1.33 


^lll, IV 


95.2 


93.7 


1.02 


80.5 


56.6 


1.42 


IV 


90.5 


76.8 


1.18 


68.3 


41.5 


1.65 


T stage 


>T3, T4 


57.1 


23.4 


2.44 


43.9 


28.8 


1.52 


Nodal metastases 


Yes 


95.2 


88.4 


1.08 


78.0 


45.3 


1.72 


N stage 


>N2b 


76.2 


49.5 


1.54 


58.5 


30.2 


1.94 


Differentiation 


Moderate, poor 


95.2 


98.9 


0.96 


100.0 


98.1 


1.02 


Poor 


66.7 


74.7 


0.89 


53.7 


34.0 


1.58 


EGFR 


Positive 


66.7 


52.6 


1.27 


78.0 


75.5 


1.03 


TILs 


Low, moderate 


90.0 


46.8 


1.92 


85.4 


84.9 


1.01 


Low 


35.0 


10.6 


3.30 


53.7 


41.5 


1.29 


Abbreviations: DR = detection rate; FPR = false-positive rate; HPV = 


= human papillomavirus; LR = likelihood ratio; OPSCC = oropharyngeal cancer; TIL = tumour-infiltrating lymphocyte. 



12.8%; compared with LR=8.6 (DR = 70.6% and FPR = 8.2%) 
using our model) and the 3-year survival was 91.5% for low-risk 
patients and 50% for high-risk patients (compared with 94% and 
36.8%, respectively, using our model). Clearly, both models have 
significant predictive power. Interestingly, we found no association 
between TIL levels and nodal status in HPV-positive OPSCC 
(P = 0.43). 

Retrospective studies have intrinsic limitations. In this study, 
cases were excluded if there was insufficient archival material for 
full pathological and immunochemical analysis. However, demo- 
graphics and survival of excluded patients did not differ 
significantly from those included in the study (Supplementary 
Information). Most patients in this study came from the same 
geographical area of England, and we cannot therefore assess the 



effect of geographical, genetic, or ethnic influences on our 
findings. The number of patients with HPV-positive/TIL low 
tumours was relatively small and, thus, any conclusions drawn 
must perhaps be treated with some caution. However, this may 
simply reflect the fact that the vast majority of HPV-positive 
tumours are associated with a prominent immune response. 
Similarly, it could be argued that the heterogenous treatments 
received by patients in this study complicates the analysis. 
However, the uniformly improved survival of HPV-positive 
patients, regardless of treatment modality, argues that the reason 
for this improved prognosis is related to an intrinsic feature of 
these tumours. Prospective analysis of our prognostic model will 
be informative, and a multicentre validation study is currently 
under development. 



498 



www.bjcancer.com | DOI:1 0.1 038/bjc.201 3.639 



Improving outcome prediction in OPSCC 
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False-positive rate (1 -specificity) False-positive rate (1 -specificity) 

Figure 4. Receiver-operating characteristic (ROC) curves for 3-year 
HPV-positive OPSCC mortality. (A) ROC curve for 3-year HPV-positive 
OPSCC mortality derived from the UHS patient cohort (training set) 
using a model containing TIL levels, smoking, and T stage. The line of 
no effect is also shown (dashed; area under the curve (AUROC, 50%). 
This prognostic model was highly predictive of HPV-positive OPSCC 
mortality, with an AUROC of 0 • 87. A cut-off score of - 0.945 was 
selected (circled in the figure) based on this ROC curve (DR, 72.7%; FPR 
,10.2%) and applied to the PFT/BLT (test) cohort. (B) ROC curve for 3- 
year HPV-positive OPSCC mortality derived from the PFT/BLT patient 
cohort (test set) using the predictive model containing TIL levels, 
smoking, and T stage derived from the UHS training set. The line of no 
effect is also shown (dashed; (AUROC, 50%)). The cut-off score of 0.945 
identified in the training set is circled in the figure. This score was highly 
predictive of 3-year mortality (AUROC, 0 • 82), with 66.7% of dead 
patients testing positive (DR = 67%) and 5.6% of alive patients testing 
positive (FPR= 1 1%). The resultant LR for prediction of 3-year mortality 
was 1 1 .9. 



In summary, we show that TIL levels predict for survival in 
OPSCC patients. High TIL levels are significantly associated with 
HPV status, suggesting that the reason for improved survival in 
most HPV-positive OPSCC is the presence of an adaptive host 
anti-tumour immune response. Survival in patients with HPV 
positive, TIL low tumours, is not significantly different than in those 
with HPV-negative disease. A prognostic model based on low TIL 
levels, heavy smoking, and late T stage is extremely effective at 
identifying a group of HPV-positive patients with poor survival. 
Tumour-infiltrating lymphocyte evaluation can be performed 
quickly on diagnostic biopsy, and more complex analysis of the 
T-cell infiltrate does not improve on the predictive power of the 
simpler scoring method. This is an effective tool for identifying 
HPV-positive OPSCC patients who are likely to respond poorly to 
treatment, and is highly relevant for the clinical and pathological 
evaluation of HPV-positive OPSCC. 
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